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Abstract 

Linearized four-derivative gravity with a general gauge fixing term is con- 
sidered. By a Legendre transform and a suitable diagonalization procedure it 
is cast into a second-order equivalent form where the nature of the physical de- 
grees of freedom, the gauge ghosts, the Weyl ghosts, and the intriguing "third 
ghosts", characteristic to higher-derivative theories, is made explicit. The sym- 
metries of the theory and the structure of the compensating Faddeev-Popov 
ghost sector exhibit non-trivial peculiarities. 
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Introduction 



Theories of gravity with terms of any order in curvatures arise as part of the low 
energy effective theories of the strings [1] and from the dynamics of quantum fields 
in a curved spacetime background [2]. Theories of second order (four-derivative the- 
ories in the following) have been studied more closely in the literature because they 
are renormalizable [3] in four dimensions. This property spurred Renormalization- 
Group studies [4-8] , including attempts to get rid of the Weyl ghosts (also known as 
"poltergeists") usually occurring in higher- derivative (HD) theories. On the practical 
side, HD gravity greatly affects the effective potential and phase transitions of scalar 
fields in curved space-time, with a wealth of astrophysical and cosmological proper- 
ties [9]. These phenomenological applications contributed to keep alive the theoretical 
interest, as illustrated by the most comprehensive introductory study available [10]. 

Besides the renormalization properties [3], all that was known about the struc- 
ture of the (classical) theory was the particle contents, as read out of the linear 
decomposition of the HD propagator into pieces with second order (particle) poles. 
Some related aspects of the equations of motion [11] were also elucidated. Definite 
theoretical progress came from a procedure, based on the Legendre transformation, 
devised to recast four-derivative gravity as an equivalent theory of second differential 
order [12]. A suitable diagonalization of the resulting theory was found later [13] that 
yields the explicit independent fields for the degrees of freedom (DOF) involved (usu- 
ally including massive Weyl ghosts), thus completing the order-reducing procedure. 
One should notice that theories with terms of higher order in curvatures have the 
same DOF and propagators of the four-derivative one, since the higher terms do no 
contribute to the linearized theory. 

An alternative order-reducing method has been proposed [14] that introduces 
an auxiliary field coupled to the Einstein tensor G^ u (or to the scalar curvature R) 
and featuring a squared term. It can be shown that this method is equivalent to the 
one based on a Legendre transformation with respect to G^ u (or to R), the auxiliary 
field being a redefinition of the "momentum" conjugate to them. 

The studies [12-14] above were carried out for the (non quantizable) Diff-invariant 
theory. An exploration of the method in the presence of gauge fixing terms has been 
done in a simplified HD gauge field model [15]. In this paper we implement this 
procedure for four- derivative gravity. 

Amongst a crowd of positive and negative norm, gauge-independent and gauge 
ghost, masless and massive states, the famous "third ghosts" arise. These subtle 
ghosts, missed in [4] and properly accounted for in [5] and there since, first emerged 
from a functional determinant in the context of Path Integral quantization. Now they 
appear as the poltergeist partners of the usual gauge ghosts. 

In Section 1 we present our starting Diff-invariant four- derivative theory. A 
very general gauge fixing term is then introduced that includes the most used ones 
as particular cases. Being interested in the propagators and in the ensuing DOF 
identification, we focus mainly on the free part of the Lagrangian. Self-interactions 
and interactions with other matter fields are embodied in a source term and may be 
treated perturbatively. Then the relevant total gauge fixed linearized theory is worked 
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out. Section 2 deals with the order- reducing procedure that leads to the diagonal- 
ized second- derivative equivalent theory. The structure of the propagators and the 
identification of the DOF is then worked out in Section 3. The Faddeev-Popov (FP) 
compensating Lagrangian is studied in Section 4, where an order-reduction of the 
fermion sector is also carried out. Particular attention is paid here to the identifica- 
tion of poles and to the striking cancellation mechanism of ghost loop contributions. 
Related to this, a discussion of the BRST symmetries involved is also made. The 
above results are summarized and further elucidated in the conclusion. 

The definitions of the spin projectors and related formulae, the basis of local 
differential operators, and the notations and conventions used throughout the paper 
have been collected in Appendix I in order to render the work almost self-contained 
and more readable. Secondary calculations regarding the conditions from locality 
on the gauge-fixing parameters and the order-reduction of the HD fermionic FP La- 
grangian have been respectively moved to two Appendices. 

1 The Linearized Lagrangian 

We consider a general four-derivative theory of gravity 

£-HD — £"inv + £g + C m , (1) 

where C m is the coupling with matter, 

= V9[aR + bR 2 + cR^R^] , (2) 

is the most general Diff-invariant gravitational Lagrangian of second order in curva- 
tures (the squared Riemann tensor has not been considered as long as a topologically 
trivial 4D spacetime is assumed so that the Gauss-Bonnet identity holds), and 

C g = ^g\xnh]G»»X u [h], (3) 

where 

X ^[h] = A^-XD^h , (4) 

Gvw = tiD p D p g^ u - i 2 -D ip D u) + + ^R^ + &Rg pu , (5) 

is a general gauge fixing term that depends on six gauge parameters and generally 
contains HD as well as lower-derivative (LD) terms. One may obtain the gauge fixings 
used in [4]- [8] by specializing these parameters. 

In order to study the propagating DOF of the theory we work the quadratic 
terms in h^ u out of Chd- Dropping total derivatives, they write 

C HD = £% v + £W + £ s = ^(P™ a + P^)h^ + jC s . (6) 

The source term C s includes the interactions with matter fields <fi and all the self- 
interactions of hfj, u affected by the Newton constant Gn- Here and in the following 
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the indices are rised and lowered by r]^ and usually omitted for simplicity whenever 
no ambiguity arises. 

The differential operator kernel for the diff-invariant part is 



. p(2) _ p(S) 



+ 66D 2 P (S) + cD 2 



:p(2) + 2 p( s ) 



The gauge fixing contribution 



yields 



4 2) = \( A " - Xd^h^a^ - i 2 d^d v + ^)(A V - \d»h) 



pa = _ nA 2 [(^-e 2 )n + e3](3P (5) + P (H/) + P {5H/} ) 

+ £ 2 D 2 pW - □ + 6] Qp (1) + 

+ An [(& - + 6] (2P (H/) + • 



(7) 



(8) 



(9) 



One recognizes in (8) the linearized an d the h- independent part of Q^ v , which 

we call in the following. 

Thus the complete HD differential operator kernel is 



pinv _|_ pg 

I n ( cn + a)p( 2 >-±n(&n + &)P (1) 



+ □ 



a + 2(36 + c)U - 3 A 2 ((& - + 6)] 

-(A-i) 2 n((6-6)n + e3)P w 
-A(A-f)n((ei-6)n + e3)P {5H/} . 



(10) 



By decomposing P in its HD and LD parts, namely 

p = Mn 2 + Nn , 



(11) 



where 



M = 



c _ P m _ I 6P (i) 

2 2 U 

+ (2(36 + c) - 3A 2 (6 - 6)) ^ (5) - (A - 1) 2 (6 - 
_ A (A-1)(6-6)^ {5W} 



(12) 



AT 



2 2 

- (a + 3A 2 6) P {S) - (A - 1) 2 £ 3 P (H ° - A(A - l)£ 3 P {SH/} 
equation (6) may be written as 



C HD = -hU(MU + N)h + £ s 



(13) 
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Dropping total derivatives, it can be given the more convenient form 



1 1 

£ HD [h, ah] = -(Uh)M{Uh) + -hN(Dh) + C s 

The HD Euler's equation takes also the form 

D(MD +NY u ' pa h pa = T» v , 

where = —4f^- . 

Otlu v 



(14) 



(15) 



2 Second order equivalent theory 

In order to perform a Lorentz-covariant Legendre transformation [13] [16] of our HD 
Lagrangian, the form of (14) trivially suggests defining the conjugate variable 



7T 



[IV 



dC 



HD 



<9(QV) 



One finds 



7T = M(Uh) + -Nh + 0(G N ) 



(16) 



(17) 



where the contributions from the gravitational interactions may be accounted for 
perturbatively in , or may be simply ignored for the analysis of the propagating 
DOF. 

As required, (17) is invertible and gives 



Oh = M- 1 



7T--Nh 

2 



= F[h,n] 



(18) 



Notice that the operators M and iV are invertible as long as gauge fixing terms have 
been introduced. Otherwise they would project into the spin-state subspace 2 © S, 
then being singular. 

The Lorentz-covariant Hamiltonian-like function is then 



H[h,ir] = irF[h,ir}- £ HD [h,F[h,ir}} 

= - \-Nh - ttI M" 1 \-Nh - 
2 [2 J [2 

The equations of motion turn out to be the system of canonical-like equations 

m 



ah = 

□tt = 



dn 

on 
~dh 



(19) 

(20) 
(21) 



The familiar negative sign one woud expect in (21) is absent because the defini- 
tion (16) involves second derivatives of the field h instead of the usual velocity [15]. 
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They may also be derived by a variational principle from the so called (now second- 
derivative) Helmholtz Lagrangian 



C H [h,n] =nnh-H[h,7i} 



(22) 



In fact from (22) one sees that (20) is the Euler's equation for 7r and (21) is the one 
for h. From (20) (which is nothing but equation (18)) one may work out ix as given 
by (17). Substituting it in (21) one recovers (15), namely the original HD equation 
of motion. 

Mixed 7r — h terms occur in (22). The diagonalization can be performed by 
defining new tilde fields such that 



h = h + 7r 

N (h ~\ 
7r = -—{h — Tr). 



(23) 



or conversely 



h = N' 1 
N- 1 



71 



Nh + n 



Nh-n 



Then Ch finally becomes the desired LD theory 

C LD = -hNn~h--7}(Nn + NM- 1 N)n + C s 
2 2 



(24) 



(25) 



where 



NM~ l N 



^_ p{2) _ Jf_ p(l) 

2c 2& 

a 2 (6-6)-3A 2 e 3 2 2(35 + c) p(sl 



+ 



2(36 + c)(6-6) 

_ ~ !) 2 ^3 p(W) 

6 -6 

_ - 1 )^3 p\SW} 

6-6 



(26) 



For further discussion, a most enlightening expression for (25) is obtained by sepa- 
rating the gauge-dependent parts 



1 



\2 J 2c 2(36 + c) 



7T 



(27) 
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where 



Gafi = + 6^a/3 = (28) 

^ - 6 6Q fa/9 + 6 > ( 29 ) 

and the form of x has been displayed in (8). 

The physical meaning is now apparent: h and tt describe the massless and the 
massive DOF of the theory respectively Notice that the gauge- invariant part of the 
kinetic term of n reproduces that of the Fierz-Pauli theory [17]. 

The LD Lagrangian (27) thus obtained is non-local for arbitrary gauge pa- 
rameters. However, we can have locality for a particular choice of parameters (see 
Appendix II). Even for this choice, an unpleasant feature of (27) is that the scalar 
subspaces S and W appear mixed as long as the transfer operator pi sw } occurs in 
N and NM^N. 



3 Linearized theory and propagators 

In order to avoid unessential complications due to the S-W mixing that obscures the 
identification of the propagating DOF, one may redefine the field h^ u as 

V = (Q~X° V (so) 

where 

Q(A) = p(2) + p(i) + ^ P W_2(A_l| p{w} (31) 

is invertible and becomes a numerical matrix for A = — 2 , namely Q(— 2) — fj — ^fj . 
This choice is compulsory if we wish to avoid polluting the source term with non- 
locality. The operator P transforms to 



P = QPQ = MD 2 + NO , (32) 



where 



M = £p(2) _ I^p(i) 



2 

2 



+ ^^5^2(36 + C)P(»'» - ±<^A(6 - 6)P' S » , (33) 

do not contain the operator pi sw } anylonger. Then equation (13) may be written as 

C HD = l -hO(MO + N)h + C s , (34) 
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or, dropping total derivatives, 

C% ) D [h,Uh] = \{Uh)M{Uh) + \hN(Uh)+fh . (35) 



The particle interpretation of (35) is now the central question. On one hand 
we can start from the HD theory (34) and, after inverting the projectors, obtain the 
quartic propagator 

AHD[h] = ( C n + a )D ^ + T (A - 1) 2 [2(36 + c)D - ap^ (36) 

2 .p(i) _ HZ p(s) 



On the other hand, the quadratic propagators arising from the LD theory (analogous 
of (25)) for the new hat fields are 

a ld [K] = A P (2)_gZ A2 P (w) 

11 an 4 (A - \) 2 aU 

2 -P {1) - 2 4t^4i^P {S) , (37) 



4(A-i) 4 6^ 

A ^r#i = 2c P (2) | 27 A2 2 ( 3fe + c ) pCHQ 

M a(cD + a) 4 (A - l) 2 a[2(36 + c)D - a] 

+ 2 6 pd) + 27 ^ (^jg) p(5) 

^ 6 + 6) 4 (A - 1) 4 ^ [(6 - + &] 

As expected, the LD quadratic propagators sum up to give the HD quartic one, 
namely 

A HD [h] = A LD [l] + A LD [ir] (38) 

Notice that if we had not performed the Q-transformation, the propagators would 
have been 

A HD [h] = QA HD [h]Q 

A LD [h] = QA LD [h]Q (39) 
A LD [n] = QA LD [ir]Q , 

with the mixing pi sw } occurring in them. 

The DOF counting can be readily done on (37). Since we are dealing with a 
properly gauge fixed four-derivative theory, all the fields in h^ v do propagate and 
therefore we have a total of 20 DOF (10 massless and 10 massive). According to the 
dimensionality of the respective spin subspaces, they are distributed as 5, 3, 1 and 1 
DOF for the spin-2, 1, 0$ and 0w respectively, which sum up 10 DOF for the massless 
fields, and the same for the massive ones. In the (massless) /i-sector, the spin-2 space 
contains the two DOF of the graviton plus three gauge DOF, and the remaining five 
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DOF are also gauge ones. The 7r-sector describes the five DOF of a spin-2 poltergeist 
with (squared) mass a/c [13] [17], one physical scalar DOF with mass — a/2(3b + c), 
three third ghosts with gauge-dependent masses £3/6 and one third ghost with mass 

6/(6-6)- 

In the absence of gauge fixing only the projectors P^ and P^ are involved, and 
pmv j n (j-j can k e inverted only in the restricted spin subspace 2 © S. In principle, in 
that case one is left with eight DOF in the LD theory, namely, the massless graviton, 
the massive spin-2 poltergeist and the physical scalar. This is generally so as long 
as critical relationships between a, b and c are avoided [10] [13] that make the order- 
reducing procedure singular. In those cases some DOF may collapse and the theory 
may have fewer than eight DOF and/or larger symmetries. Fast DOF-counting recipes 
for gauge theories can be found in [18]. In ordinary two-derivative gravity, each of 
the four gauge-group local parameters of Diff-invariance accounts for the killing of 
two DOF, leaving the two DOF of the graviton out of the ten DOF of . In four- 
derivative gravity each gauge-group parameter instead kills three DOF so that the 
initial twenty DOF reduce to the eight DOF quoted above. The mechanism is well 
illustrated in four-derivative QED [15], where one initially has eight DOF and the 
gauge invariance suppress three of them, leaving one photon and one massive spin-1 
poltergeist. 

As detailed in Appendix II, the free part of the Q-transformed LD theory can be 
made local for a particular choice of gauge parameters. However, using (5(A), for A = 
6/(46 +c), moves the non-locality to the source term, namely to the interactions. This 
can also be avoided by requiring that A = — 2 , in which case Q becomes a numerical 
matrix, but this gives rise to a condition on the parameters b and c of the starting 
Diff-invariant theory. Leaving aside the interpretation of such a restriction, this does 
not mean that we had anyway obtained a sum of independent Lagrangian theories for 
each (massless and massive, spin 2, 1, O5 and 0w) particle, notwithstanding the fact 
that the spin subspaces appear well separated. This is not posible, as illustrated for 
instance by the fact [19] that there is no second-order tensor local theory for spin-1 
fields. 



4 Faddeev-Popov compensating terms 



As usual, the gauge fixing term (3) together with the compensating (HD) Faddeev- 
Popov Lagrangian can be expressed as a coboundary in the BRST cohomology, 
namely 

(40) 



Cg + CpP 



where C are FP fermion ghosts and B is an auxiliary commuting field. 
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In order to study the propagators of the new fields C and £>, it suffices to 
consider the linearized objects 

Gap = fl^o/jD ~ ^dadp + fgTfep (41) 

= (6n + 6)^ + [(6-6)n + 6]^a/3 , 

and the BRST symmetry given by the (linearized) Slavnov transformations 



qIi — T) C a 

sC - a = ° (42) 
sC a = B a y ' 

sB a = , 

where 

£W = Vtfdv + Vitpdn (43) 
is the gauge symmetry generator. With the diagonalization 

B a = B a - X a [h] (44) 

the linearized (40) becomes 

4 2) + C fF = \x a [h]G a h ,V[h] + C a g ayX ^V^C? - \B-Q^ , (45) 



ah — T) r 10 



and one has 

sC a = 

sC a = B a - X a [h] 1 } 

sB a = x^V^C? . 

Of course, they reflect the trivially Abelian gauge symmetry group G to which the 
infinitesimal diffeomorphisms reduce in the linearized theory. In the complete non- 
polynomial theory there are couplings between h^ v , C , C and B , and the non- 
Abelian symmetry yields a more complicated set of ^-transformations in which, for 
instance sC a ^ 0. One should also notice that the Q-transformation used in Section 
3 is ininfluent in the form of Cfp, as can be checked by computing it with the 
corresponding operators V = Q~ l V and x a>iv — (xQ) aflI/ - 

The fermionic sector of the FP Lagrangian above, namely 

C HD [CC] = C a [(&□ + Z 3 )ne af) + 2(1 - A)((6 - + 6) □"ad C p (47) 

is higher-derivative whereas, in constrast with ordinary second-order theories, now 
the auxiliary bosonic field B does propagate according to the Lagrangian 

C[b] = ~B a + 6)^ + [(6 - + &] ^] ^ (48) 
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which is already LD and always local. We can also perform an order-reduction of 
C HD [CC] (Appendix III), yielding 

C LD [EEFF] = E a ^ 3 9 a/3 + 2(l-\)^u; a/3 )nE 13 (49) 

- F" (kfan + ^) 9afi + 2( *~^ 3 ((6 - + fc)^) F' 3 , 

where E a + F a = C a and E a + F a = C a . The Lagrangian (49) is local for the same 
choice (see equations (102)- (104) in Appendix II) of gauge parameters that makes 
Cld in (25) local. 

From (48) one directly reads 

^ = 55T& + (&-&)□ + 6 ' (50) 

whereas the HD (oriented) propagator 

A HD^ n] , u 1 / m n 

[UUJ " &□ + &)□ 2(1 -A) ((&-&)□ + &)□ 1 J 



6 \u Zin + &J 2(i-A)6\n (&-&)□ + &/ 

obtained from (47) splits into the (also oriented) LD propagators 

- w + (52) 

AlPF] = "&(&□ + &)"" 2(1 -A)&«6 -&)□+«" ' <53) 

that can be derived from (49) as well. 

In (52) one counts four FP fermion ghosts E and four .E with massless poles, 
giving eight negative loop contibutions that compensate for the eight massless gauge 
ghosts quoted in Section 3. The compensation of the third ghosts contains non trivial 
features which are characteristic to the HD theories. From (53) one has that the FP 
fermion ghosts F and F give six negative loop contributions with propagator poles 
at £3/^1 and two at £3/(^1 — £2)- This over-compensates the (three plus one) third 
ghosts. Here is where the new boson FP ghosts B, propagating with (50), come to 
the rescue: they give three positive contributions with poles at £3/6 and one at 
£3/(6 _ £2), thus providing the complete cancellation of ghost loop contributions. 

This matching of the ghost masses is a consequence of the interplay of the 
order-reducing and BRST procedures. The master relationship is 

gW- 1 = gih)- 1 _ gfr)- 1 ^ ( 54 ) 

where the massive poles of the third ghosts are displayed by 1 and 1 , the 
latter also having massless zero-modes. To see it we find useful defining the differential 
operator 

Z% = x^,^ = □ + 2(1 - A)a£l , (55) 
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and the differential kernels 

K^^G^Z} (i = M,7f) (56) 

occurring in the FP Lagrangians above and worked out in (47) and (49). The poles 

of the (massless) gauge ghosts lie in Z~ l whereas the operator has no poles. 

From (54) and (56) it follows that 

KW- 1 = K&)- 1 _ K&r 1 (57) 

which also reads A HD [CC] = A[EE] + A[FF] . Thus the ^-fields and the F-fields 
inherit the massless and the massive poles respectively. On the other hand A [B] = 
—Q^ 1 , so that also the boson B-fields share the same massive poles. 



The symmetries of C^ D + Cp^> are not trivial either. The symmetry of the 



'LD 



(invariant part of the) HD theory under the group G of the gauge variations 



&hp V Dpv,a£ 



(58) 



is inherited by the LD theory via (17) and (24) with the variations 



Sh 



pV 



Sfpu = 



IpU 



aiN-'MYZ 



[IV 

£l D p(l)pa , (6 ~ 6) n p(W)pa 

( I'" " t " 

?3 43 



+ 

V 



^ p<r,a° 



£l D + 6 p(l)pa , (6 - 6) D + & p(W)pu 
« <53 



pv 



(59) 



pa,a L 



pv 



1 - y pa,a^ 



both depending on the same four gauge-group parameters e a (x). One may check that 
Sh^ + 6^^ = Sh^. However, the free invariant part of the LD theory (25) exhibits a 
larger symmetry group, namely the fields h and tt may be given independent variations 



pv *D pv,a£ 



la 
l/a 



(60) 
(61) 



thus doubling the number of group parameters, with the original symmetry as a 
diagonal-like subgroup G\ C G x G , which is isomorphic to G [15]. One may then 
look at 

Ci D [h] = \x[h}Q {K) x[h] (62) 



and 



1 



(63) 



occurring in (27), as separate gauge fixings for the symmetries (60) and (61) respec- 
tively, and wonder what happens with the whole BRST scheme. 
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The separate S'-transformations would be 



sE a = 



V E a 



S7l„ 



1 fiu 

sF a = 



V F a 



sE a = B' a - x a [h] 

sB ,a = x ^V^pEP 

so we are led to write 



FP 



— S 



sF a = 
sB" a 

1 



B 



Ha 



X a [n] 



(64) 



(65) 



+ s 

1 

2 

+ E a gj$ x p ^v^i 



x a [h}g ( Sx"[h] - -\ 



1 

2 

1 

2' 

£0 



F <*gW x ^v^Ft 



(66) 



_ + \b"^b" 



Thus (49) agrees with the fermionic sector of (66). 

Equations (64) define twocohomologies {V; s} and {V; s} with cohomological 
spaces V = {h, E, E, B'} and V = {tt, F, F, B"} respectively, both being copies of 
the original V = {h, C, C, B} of (46) with boundary operator s. The polynomial (65) 
is then an exact cochain in the cohomology {V; s} © {V; s} = {V © V; s © s}. 

The cohomology characterizing the HD theory appears as a subcohomology 
{Vi; si} of the direct sum above. The subspace V± C V © V is defined by 



B' a = OfB 



F a = OfC® 



per 



B" a 



= O'^B 13 



where 



(67) 

(68) 
(69) 

-a/; (70) 

6 (71) 

are invertible linear operators satisfying O' + O" = 5, and Si is the restriction to V\ 
of s © s. Then this subcohomology is nothing but the original one {V; s} of the HD 
theory since (67) defines an isomorfism V A V\ and s\ becomes s , that is 



/nlla 

U I3 



-u(n- 1 my; u 

—00 H 




sE a + sF a 
sE a + sF a 
s5 to + s5" Q 



= sC a 
= sC a 
= sB a 



(72) 
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as a consequence of (67) and (68)- (71). In other words we have i 1 1 o s © s o i 1 = s. 
Moreover, we recover the Lagrangian (48) for B, namely 

C*[B'B"\ = --B^Q^B'P + -B" a g$jB" fi 

= -\B a O'2Q { ^O^B p + l -B a O'^g^Ol P B p (73) 

= —B-g^ = C[B\ 

The subgroup G\ C G x G, associated to {T4; si} and isomorphic to G, is obtained 
by taking the group parameters e' and e" as functions of four parameters e by means 
of the equations 

V^ a e' a = 0%V pa , a s a (74) 

T) r" a n" pa T> r a 

These are derived by imposing the relations stemming from (59) on the otherwise 
independent variations (60) and (61), and yield 

e" a = Ofs? , (76) 

so that e = e' + e" . The subgroup G\ is, by definition, a symmetry of the whole (non 
gauge-fixed) Lagrangian (and also separately of the interaction terms) since one has 
5h^ u + Stt^i/ = 5h^ u , whereas G x G is broken by the interaction terms. 

5 Conclusion 

The interplay of gauge invariance and higher differential order in field theories gives 
rise to a remarkable diversity of particle-like states which are encripted in the orig- 
inal field variables. In the four- derivative tensor theory of gravity here studied, the 
doubling of the initial conditions for the (fourth differential-order) equations of mo- 
tion translates into a doubling of the effective number of particle-like DOF obeying 
second differential-order evolution equations. They describe physical (positive Fock- 
space norm) states together with an outburst of massless and massive ghostly states 
which are unphysical because of their negative norm and/or gauge dependence. Be- 
yond the methodological interest, its analysis provides an enlarged context for the 
traditional gauge theories and BRST symmetries of physical relevance, which also 
enlightens the nature of some states already encountered in former classical works on 
higher- derivative gravity. 

Four-derivative gravity is particularly interesting to study as long as the empha- 
sis traditionally given to its applications has overlooked many details of its theorical 
structure. Amongst the particle-like states of the gauge-fixed theory, there are phys- 
ical ones (a massless graviton and one scalar, reminiscent of the Brans-Dicke field), a 
massive spin-2 gauge independent Weyl ghost (unphysical norm), and two families of 
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gauge-dependent fields: the usual massless gauge ghosts and the novel massive third 
ghosts. This elusive new breed of ghosts firstly arose in the exponentiation of the 
functional determinant of the differential operator Q^ h \ 

In the presence of (generally HD) gauge fixing terms and of the correspond- 
ing compensating FP Lagrangian, the order-reducing procedure reveals remarkable 
features of the underlying BRST symmetry associated to the four-parameter gauge 
group G of infinitesimal diffeomorphisms. In parallel with the doubling of the fields, 
there is a doubling of the gauge symmetry of the free part of the (second-order) LD 
equivalent theory. Out of this G xG larger symmetry, both the interaction terms and 
the consistency of the BRST algebra select a diagonal-like subgroup Gi , isomorphic 
to G , as the only symmetry of the complete LD theory, in agreement with the ocur- 
rence of Diff-invariance as the only symmetry of the starting HD theory. However, 
restricting ourselves to the free LD theory and considering its G x G symmetry, the 
LD gauge-dependent terms can be viewed as separate gauge fixings for both group 
factors. The (gauge-dependent) unphysical propagating DOF so introduced then ap- 
pear as the respective gauge ghosts, which are massless for the first group factor and 
massive for the second one, thus giving further meaning to the famous third ghosts. 
Moreover, the separate symmetry of the the gauge-independent part of the physical 
and poltergeist sectors of the LD theory illustrates also how their kinetic terms repro- 
duce the structure of the Einstein's and Fierz-Pauli theory, thus describing (massless 
and massive) spin-2 fields respectively. 

In correspondence with the appearance of a new class of massive gauge ghosts, 
the compensating FP Lagrangian also contains a greater number of propagating fiels. 
These come from the HD doubling of the FP anticommuting fermion fields and from 
the boson fields, which are just auxiliary decoupled artifacts in ordinary two-derivative 
theories and now propagate and couple to the gauge-independent fields. The negative 
loop contributions of the massive fermion FP fields yield twice the amount needed 
to compensate for the third ghost loops, and it is just the positive contributions 
of the boson FP fields that provide the exact balance. This striking compensation 
mechanism, peculiar to HD gauge theories and easy to extrapolate to higher than four- 
derivative theories, illustrates well the power and richness of the BRST procedure. 
Of course, checking the exact cancellation of ghost loop contributions would require 
considering the actual residues of the propagators and vertex couplings arising in the 
complete non-polynomial theory, a task which is beyond the purposes of this work. 

A final comment on locality is in order. From a HD local theory, the order- 
reducing procedure leads to an equivalent two-derivative theory. For scalar theories, 
the LD counterpart is directly local [16]. In gauge theories of vector fields there is 
always a choice of the gauge fixing parameters for which it is also local [15]. For 
tensor fields, the example studied in this paper tells us that obtaining an equivalent 
LD theory with independent free Lagrangians for the different spin states is not 
compatible with locality, although one comes close to this goal by suitably picking 
the gauge parameters. This obstruction is related to the more complex structure of 
the constraints of the tensor field theories, like the one that prevents having a tensor 
local theory of second differential-orden for spin-1 fields. 
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Appendix I 

We use the notations 



a» = d v w u 



h = K 

and the Minkowsky metric is rj^ = diag(-\ ). 

The spin projectors are 

P fv,po = -^li{pQvo) ~ T^iivQpo (77) 

Pj$** = (78) 

PjlV,pO = ~^\iv®po (79) 

PZl = <W (80) 
They are symmetric under the interchanges 

/i <-> v , p <-> <T , pu ^ pa , (81) 

idempotent, orthogonal to each other, and sum up to the identity operator in the 
space of symmetric two-tensors, namely 

W = 5 = (P (2) + P (1) + P (S) + P (W) U P « ■ (82) 

These projectors are constructed using the transverse and longitudinal projectors for 
vectors fields 

<«V = ^ • (84) 

We also use the transfer operators 

P$2 = <V«V (85) 

= (86) 

from which we define 

p{SW} = p{WS} _ p(SW) + p(WS) _ ^ g7 ^ 
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They have non-zero products 



We define also 



p(SW)p{WS) 
p(WS)p(SW) 
p(S) p(SW) 
p(W) p(WS) 
p{SW}p{SW} 
p(S) p{SW} 
p(W)p{SW} 



3P (S) 

3 p(W) 

p(SW)p(W) _ p(SW) 
p(WS)p(S) = p(WS) 

3 (p(S) + p W ) 

p{SW}p(W) _ p(SW) 
p{SW}p(S) = p(WS) 



I [IV, pa 



= V„V P a = 3P (5) + P {W) + P {SW} 



(89) 
(90) 
(91) 
(92) 
(93) 
(94) 

(95) 



Here we collect some formulae which are useful for dealing with combinations 
of the operators above. Inverse: 

M = \ 2 P {2) + AiPW + Xs pW + Xw P<W + \ SW P {SW} 

M -l = J_p(2) + J_ p(l) + V p (S) + X S p(W) 

A2 Ai AsAvk — 3A§ W A^A^ — 3A| W 

A^Apf — 3\ 2 SW 

When computing symmetric products like (26), operators in the subspace S © W of 
the form 

Ufa, r w , T SW ) = r s P (S) + T W PW + r sw P {SW} (96) 
occur, for which one has the product law 

Q(t~s, t~w, t~sw)Q(^s, ^w, ^sw)Q(ts, t~w, t~sw) = (97) 
Q( r s^s + 3Tg W \ w + 6t s t S w^sw , 

T w^w + 3ts W \ s + 6t w t sw \ sw , 

l~Sl~W^SW + TSTSwAs + TswTw^W + ^sW^Sw) 



A basis of zeroth differential-order local operators with the symmetries (81) is 
provided by fj and fj. Local second-order operators can be expanded in the basis 

C 2 n := (^P« + 3P^)n 

C 3 D := (P^ + 6P^)n (98) 

C 4 n := (P (H/) -3P( s ))n 
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Thus, a general local LD operator has the form 

4 

Q LD = Y, a ^i D + a ^ + 02V • (99) 
1=1 

Einstein's (linearized) theory displays the operator C\. Fierz-Pauli's has the 
same kinetic term and a mass term built with 77 — 77. When using the field basis 
obtained by Q(A)-transforming the theory, this kinetic term displays the operator 

Q(A)GDQ(A) = Qp< 2 > - ^ ( -^P-P {W) ) □ (100) 

For A = —2 it becomes local again, namely (|p( 2 ) - □ = dD - |C 3 D which 

describes (linearized) gravity as properly as the original operator C\U did. 



Appendix II 

The kinetic terms for h and tc in (25) contain the operator NO which is local for 
arbitrary gauge parameters. In fact one has that 

N — ad — £A - A(A - 1)&C 3 - 6(A - 1) 2 C 4 (101) 

However, the "mass term" — ^ttNM^Ntt is local only for a choice of gauge parame- 
ters obeying the conditions 

6 = -c| , (102) 
a!6^6 _ _36 + c 

2c CI ~ 46 + c ' UUdj 

A = — ^— , (104) 
Ab + c v 7 

that are obtained by requiring NM~ l N to be a linear combination of fj and 77. This 
leaves one of the parameters £ still arbitrary. This same conditions make the fermion 
Lagrangian (48) local. 

In view of the conditions above, a theory in which 46+ c = does not have a local 
LD equivalent. But this case was critical already for the complete Diff- invariant theory 
since it is not regular in R and a general-covariant Legendre transform (see equation 
(8) of [13]) cannot be performed. With gauge fixing terms and for the linearized field, 
the (just Lorentz-covariant) Legendre transform can always be carried out (equation 
(17) is not singular for 46 + c = 0) and we instead have non-locality of the LD theory. 
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When we consider the Q-transformed theory, the potentially non-local operator 

is 

NM-'N = -P( 2 ) + - (A ~ 1)2 , f - PW 
2c 27 A 2 2(36 + c) 

_ $ p(i) 4 (A - 1) 4 g (g) 

26 27 A 2 6-6 • ( ° 5) 

As explained before, we must take A = — 2 in order to keep the locality of the source 
term. In that case NO remains local. Requiring locality for NM~ X N leads to 

6 = -^6 (106) 
6 2 = ^6 (107) 

c=--b (108) 

These conditions are the same found before, but now (104) yields a condition, namely 
(108), on the parameters of the original gauge-invariant theory. 



Appendix III 

We briefly outline the order-reduction of the higher- derivative FP Lagrangian for 
anticommuting fields 

Chd = C, (□(<*!□ + h) 9^ + n(a 2 n + b 2 ) un C v + CC, + (109) 

where ( and ( are external sources which are also anticommuting. Dropping total 
spacetime derivatives, conjugate momenta may be defined as the left derivatives 

V» = ^MR = M ^uC u + \m^C v (110) 
u n z 

V» = ^^ = -M^nC v -\u» v C v (111) 



where 
from which 



M=a 1 6 + a 2 uj , U = b l e + b 2 uo , (112) 



□ C„ = M2(p»-\aT">C p ) (113) 



1 



uC, = -M~l V v + -zM vp C p (114) 



2 
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Then the "Hamiltonian" is 

H = (nO)V + (nC)V - C 

= -(p + ltfc^M-'fr-^tfC^-eCr-C^ (115) 

With the field redefinition 

C = E + F C = E + F 

<p = \M{E-F) V = \M(F-E) (116) 
the Helmholtz Lagrangian 

C H = (nC)V + (nC)V-n (117) 

becomes 

C LD = EUuE-f(Mu + MM- 1 M)f 

+ aE + F) + (E + F)( . (118) 
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